Scaling properties at the interface between different critical subsystems: The 

Ashkin- Teller model 
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We consider two critical semi-infinite subsystems with different critical exponents and couple them 
through their surfaces. The critical behavior at the interface, influenced by the critical fluctuations 
of the two subsystems, can be quite rich. In order to examine the various possibilities, we study 
a system composed of two coupled Ashkin- Teller models with different four-spin couplings e, on 
the two sides of the junction. By varying e, some bulk and surface critical exponents of the two 
subsystems are continuously modified, which in turn changes the interface critical behavior. In 
particular we study the marginal situation, for which magnetic critical exponents at the interface 
vary continuously with the strength of the interaction parameter. The behavior expected from 
scaling arguments is checked by density matrix renormalization group calculations. 



I. INTRODUCTION 

Realistic systems have a finite extent and, when they 
display a second-order phase transition, the critical be- 
havior in the boundary region is generally different from 
that in the bulk. The characteristic size of this region 
is given by the correlation length, which becomes diver- 
gent as the critical temperature Tc is approached. In 
the vicinity of the critical point, the singularities of local 
quantities, such as the surface magnetization, are char- 
acterized by critical exponents, which are generally dif- 
ferent from their bulk counterparts. This type of local 
critical phenomena has been thoroughly studied in the 
case of a free surface through exact, field-theoretical and 
numerical methods.^'-'^ 

When a system is in contact through its boundary with 
another system, the environment can influence the local 
critical behavior at the interface. If, however, the critical 
temperature of the environment, T^, is different from 
the nature of the transitions at the interface is expected 
to be the same as for a surface<^ If the environment has 
the higher critical temperature T'^> it stays ordered 
at Tc and the interface transition has the same proper- 
ties as the extraordinary surface transition.i^i^ In the 
opposite case, for < Tc, the environment is disordered 
at Tc and the interface transition is actually an ordinary 
surface transition. ii^i^ 

Here we consider the more complex problem when the 
two subsystems in contact have the same critical temper- 
ature but not the same set of critical exponents. Thus, 
the competition between the two different bulk and sur- 
face critical behaviors may result in a completely new 
type of interface critical phenomena. This problem has 
already been addressed in Ref. ^ in which the analyti- 
cal mean-field solution, in terms of ip*^ field theories, has 
been obtained and generalized by using phenomenologi- 
cal scaling considerations. Monte Carlo simulations have 



also been performed in two dimensions for interfaces be- 
tween subsystems belonging to the universality classes 
of the Ising model, the three-state and four-state Potts 
models. 

In all these examples the stable fixed points are related 
to surface critical behavior and the expected renormal- 
ization group (RG) phase diagram is the one given in 
the upper part of Fig. O For weak interface couplings, 
the junction renormalizes to a cut, and we have the 
same local critical behavior as for a free surface, whereas 
for strong couplings, the interface becomes ordered at 
the bulk transition temperature. For some intermediate 
value of the couplings, there is a special interface transi- 
tion fixed point, involving new critical exponents, which, 
however, can be expressed in terms of the bulk and sur- 
face exponents of the two subsystemsi^ 

In the present work our purpose is to examine the dif- 
ferent types of possible interface critical behavior which 
can be realized. Thus, we consider situations where a 
weak interface coupling can be irrelevant, relevant or even 
truly marginal. We are particularly interested in the lat- 
ter case. A convenient system, for which all these dif- 
ferent situations can be realized, is the two-dimensional 
(2D) Ashkin- Teller (AT) model, ^ or its one-dimensional 
(ID) quantum version.^'* 

By introducing two Ising variables per site, the AT 
Hamiltonian can be rewritten as two Ising Hamiltonians 
coupled through a four-spin interaction,^ which is a truly 
marginal operator. As a consequence, some bulk and 
surface critical exponents are continuously varying func- 
tions of the strength of the four-spin coupling e. These 
critical exponents are known exactly through conformal 
invariance^" and Coulomb-gas mapping. 

The composite system which we consider consists of 
two AT models with the same critical temperature but 
different four-spin couplings, and thus different sets of 
critical exponents. We couple these subsystems through 
their surface spins and study the critical properties at 
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the interface while varying the strength of the interface 
coupling. We first classify the possible interface critical 
behaviors through scaling considerations, which are then 
confronted with the results of extensive numerical calcu- 
lations using the density matrix renormalization group 
(DMRG). 

The structure of the paper is as follows. The AT model 
and its basic critical properties are described in Sec. [ill 
We define the composite system and discuss its possi- 
ble interface RG phase diagrams in Sec. IIIII Results of 
numerical calculations are presented in Sec. IIVI and dis- 
cussed in Sec. El 



II. ASHKIN-TELLER MODEL AND ITS 
CRITICAL PROPERTIES 

The AT model is defined in terms of two sets of Ising 
spin variables cr, = ±1 and Ti = ±1, attached to each 
lattice site i. The usual Ising interaction K{ai<jj -f TiT-j) 
between nearest- neighbor sites i and j is supplemented by 
a four-spin interaction K^aiajTiTj, which is parametrized 
as = Ke. This latter term represents the product 
of the energy densities in the two Ising systems. We 
consider the system on a square lattice and work with the 
row-to-row transfer matrix T. In the Hamiltonian limit, 
the transfer matrix can be written as T ^ exp(— kTY^t^), 
where k is the lattice spacing in the "time" direction and 
Ti.AT is the ID quantum Hamiltonian given by 
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h is the strength of the transverse field, which plays the 
role of the temperature in the classical system. One can 
introduce a set of dual Pauli operators ^^'^ '^i+1/2 

such that 
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When the Hamiltonian in Eq. ([T]) is rewritten in terms 
of the dual variables, the couplings Ji = 1 and the trans- 
verse fields hi = h exchange their roles. Consequently, 
the homogeneous system is self-dual and the self-duality 
line is located at /ic = 1- For —1/^/2 < e < 1 this 
is just the critical line separating the ferromagnetic and 
the paramagnetic phases of the system. In the region 
-1 < e < -I/V2, for h~{e) < h < /i+(e), there is a 
so-called critical fan in which the system stays critical.— 
At the critical point, the excitation energy AE and the 
wave vector k are linearly related, AE = v^k, and the 
sound velocity is given byi^ 



TT sin(arccos e) 
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In the critical system, the basic operators are the 
magnetization Om(i) = erf (rf), the energy density 
Oe{i) = cr-cr^+i {t^t^+i) or, through duality, erf (rf ), 
and the polarization Op{i) — ufr^ . The connected crit- 
ical correlation functions display a power-law decay, so 
that + r)) - + r)) ^ r^^-^, 

where Xa is the anomalous dimension of Oa- Similarly, 
surface-to-surface correlations involve the corresponding 
surface dimensions . 

The critical properties of the AT model are exactly 
known through conformal invariance and Coulomb-gas 
mapping. The anomalous dimensions of bulk operators 
are given byi 
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The correlation length critical exponent is related to the 
dimension of the energy density hy v = 1 / {2 ~ Xe) when 
— 1/\/2 < e < 1, whereas it is formally infinite in the crit- 
ical fan. True marginal behavior implies that the scaling 
dimension of the operator 04(1) associated with the four- 
spin interaction a? a^^iT^ t^,^-^ keeps the constant value 
X4 = 2, the same as for the two decoupled Ising chains. 

The corresponding anomalous dimensions for surface 
operators an 
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One may notice that the anomalous dimensions, which 
are e dependent in the bulk due to the presence of the 
marginal four-spin interactions, remain constant at the 
surface and vice versa. 



III. COMPOSITE SYSTEM AND 
RENORMALIZATION GROUP PHASE 
DIAGRAMS 

A. Ladder and chain junctions 

A composite AT system is obtained by coupling two 
different semi-infinite subsystems through their surface 
spins. These subsystems have the same nearest-neighbor 
coupling, thus the same critical temperature. They have 
different values of the four-spin couplings e^^^ (e^^^) for 
z < (z > 0) with e(-) < e(+). 

The junction can be of two different kinds, ladder or 
chain junctional (see Fig. 6.1 of Ref. ITsI ). In the ladder 
junction, there are nearest-neighbor as well as four-spin 
couplings between sites at i = — 1 (boundary of the z < 
subsystem) and i = 1 (boundary of the z > subsystem) . 
In the Hamiltonian limit, this corresponds to a term 

V_i,i = - J((TliaJ + rliTf + eintaliafrlirf ) , (6) 

and the complete Hamiltonian is written as: 



n = n 



AT ' ^AT 



nv4 + v_i 



(7) 
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In the case of the chain junction, we introduce an extra 
line of spins at z = 0, which are connected horizontaUy to 
the two subsystems through the respective bulk couplings 
and there is a two-spin interaction associated with the 
junction in the vertical direction. In the Hamiltonian 
limit, the different terms in are extended up to i = 
and the junction involves a transverse-field term 

Vo = ~h{(^o + < + eint^^o'') • (8) 

The ladder and chain defects are transformed into each 
other through duality. In the following, we study the 
ladder problem as defined in Eq. ([7]). 



two junctions are of the form given in Eq. ([8]) and placed 
at i = and i = L/2. 

In the cylinder geometry, the first gap of Hcyi scales 
as 1 / L for a critical system, and the prefactor is propor- 
tional to the anomalous dimension of the magnetization 
at the junctior^i^ 

E,~E,^ ^xj^* . (10) 

Other local exponents are similarly related to higher 
gaps. 

Before calculating the anomalous dimensions numer- 
ically, we first consider the possible phase diagrams by 
studying the stability of the different fixed points. 



B. Basic quantities 

1. Matrix elements 

We are interested in the local critical behavior of the 
system; in particular, we want to determine the anoma- 
lous dimensions associated with the interface, for 
the magnetization density 0^(0) and xjf' for the energy 
density Oe{Q)- These can be deduced from the finite-size 
scaling of the singular part of the corresponding matrix 
elements: 

mint(i) = (0|a-(±l)|0) 
eUL) ^ (0|a-(±l)|0)-eL^i-^"", 
eUL) = {Qy-{~l)a^{m~el,^L-<\ (9) 

For the magnetization density, symmetry-breaking 
boundary conditions are needed. |0) is the ground state 
of the Hamiltonian in Eq. ([7]) and e^^^l is the limiting 
value of the interface energy density in the infinite sys- 
tem. We note that the two first matrix elements can 
be calculated on each side of the interface and there are 
two possible definitions for the energy density, ef,j^ and 
efjjt, corresponding to vertical and horizontal bonds in 
the classical model. 



2. Gaps 

These exponents can also be obtained by using confor- 
mal invarianceJ^ The classical system composed of two 
semi-infinite planes coupled by one junction is mapped 
through the logarithmic transformation into two infinite 
strips, each with a width L/2, coupled together by two 
parallel junctions at their boundaries and thus building a 
cylinder. In the extreme anisotropic limit, a Hamiltonian 
Ticyh similar to H. in Eq. ([7]), is associated to the trans- 
fer matrix along the cylinder but with two junctions and 
periodic boundary conditions. For a ladder defect, the 
two junctions are of the form given in Eq. the first 
between sites i — —\ and z = 1 and the second between 
sites i = —L/2 and i = L/2. For a chain junction, the 



C. Two identical subsystems 

We start with the symmetrical model where £*■"■' = 
e'^+) = e. In this case there are three fixed points, lo- 
cated at J = 0, J = 1 and J = oo, and corresponding re- 
spectively to two disjoint semi-infinite systems (ordinary 
interface transition), to the homogeneous system (bulk 
transition), and to a system with an ordered interface 
(extraordinary interface transition) i^^i^^i^^ 

1. Ordinary interface fixed point 

At the ordinary interface fixed point the perturbation 
takes the form JO™(-l)0™(l)-h Jei„tOp(-l)Op(l). The 
first operator, involving the product of two surface mag- 
netization operators, has the dimension 

Xint — ^ + x[^^ = 2a;^ , (11) 

and, thus, the scaling exponent of J is 

V\nt = rfint - 2:int = 1 - 2a;,^„ , (12) 

where dint = d — 1 is the dimension of the interface. This 
type of perturbation is irrelevant for ?/i„t < 0, i.e., for 
a;^ > 1/2, which happens for e > 0, whereas it is relevant 
for e < 0. The marginality condition is satisfied for e = 0, 
which is the Ising limit. The second operator, containing 
the product of two surface polarization operators, has the 
dimension jifint — 2Xp = 2; therefore, this perturbation is 
always irrelevant. 

2. Bulk fixed point 

The perturbation to the bulk fixed point introduced 
by the junction now takes the form AC'e(— l)-t- AC'4(— 1) 
where A = J — 1 and A — Je-mt — £■ The dimension of the 
first operator is Xint — Xe', thus, the scaling dimension of 
A is 

Vint = dint - Xe = d - 1 ~ Xe ^ ly^^ - I . (13) 



4 





8>0 













1 


oo 





I 


oo 













— 




— 0— — 


^ J 









01 


B 


EI 


01 


SI 


EI 



£<0 

1 oo 

o > • < o J 

01 B EI 

FIG. 1: Schematic RG phase diagram at a ladder defect with 
coupling J in the critical AT model. The RG flow is different 
for different signs of the bulk four-spin coupling e. When 
e > 0, the bulk fixed point (B) is unstable; the flow is toward 
the ordinary interface fixed point (OI) when J < 1 and the 
extraordinary interface fixed point (EI) when J > 1. When 
e < 0, the flow is reversed and the bulk fixed point is always 
stable for < J < oo. 

This perturbation is relevant (irrelevant) for v < 1 
{ly > 1), i.e., for e > (e < 0). The marginal situa- 
tion corresponds once more to the Ising limit e = 0. The 
second operator ©4 has the scaling dimension X4 = 2. 
It follows that A has the scaling dimension yint = — 1. 
Thus, the four-spin interface perturbation is always ir- 
relevant as for the ordinary interface fixed point. 



3. Extraordinary interface fixed point 

The stability of this fixed point is related to that of 
the ordinary interface fixed point. Let us consider the 
chain junction in Eq. ([5]). The ordered interface can be 
realized by setting the transverse field at the fixed point 
value h = 0. Under the duality transformation in Eq. ([2]), 
the (weak) chain junction is transformed into a (weak) 
ladder junction; consequently, to decide about the stabil- 
ity of the corresponding fixed point, one can repeat the 
argument of Sec. IIII C II 



4-. Renormalization group phase diagram 

Based on the stability analysis of the fixed points, 
the expected interface RG phase diagram is depicted in 

Fig.m 

When e < and for any interface coupling < J < cxd, 
the behavior at the interface is expected to be governed 
by the bulk fixed point. Then, the first gap in the spec- 
trum of the conformal Hamiltonian Ticyi has a 1/L de- 
pendence with a prefactor which, according to the gap- 
exponent relation in Eq. (|10p . is proportional to Xm- 



e'"'+e"<0 

J, 00 

o *•* o J 

01 SI EI 

FIG. 2: Schematic RG phase diagram at the interface between 
two different critical semi-infinite AT models with four-spin 
couplings e'~' < The interface coupling J is of the 

ladder type. The direction of the RG flow depends on the 
sign of e'"' When e'~' -|-e'^' > 0, the special interface 

fixed point (SI) is unstable; the flow is toward the ordinary 
interface fixed point (OI) when J < Js and the extraordinary 
interface fixed point (EI) when J > Js- When e^~^ -t-e^"*"^ < 0, 
the fiow is reversed and the special interface fixed point is 
always stable for < J < 00. 

On the contrary, for e > the bulk fixed point is un- 
stable. For weak couplings, J < 1, the interface renor- 
malizes to a cut and the critical behavior is the same 
as at a free surface. The first gap in the spectrum of 
Ticyi for small J can be estimated perturbatively, as in 
Sec. IIII C II It behaves as the product of the two surface 
magnetizations, vanishing as ~ L^^^™ , which is faster 
than 1/L since xf^ > 1/2 according to Eq. ([5]). This in- 
dicates that the system is asymptotically breaking into 
two pieces. For strong couplings J > 1, the interface re- 
mains ordered at the critical temperature and, through 
duality, the gap has also the size dependence ~ L"^^™, 
which is faster than 1/L. It corresponds to a vanishing 
amplitude in Eq. pop and, thus, to a vanishing interface 
magnetic exponent, a value which is linked with the local 
order at the critical point. 

In the limit e = 0, i.e., when the AT model becomes 
a system of two noninteracting Ising models, the inter- 
face coupling J is a marginal perturbation and the local 
magnetization exponent is J dependent 

a;™ (^) = 4arctan2(l/J), 4"* = 1, e = 0. (14) 

Similarly, for a chain junction,^** the local magnetization 
exponent is h dependent, 

x^rW = A arctan2(/J), x'"* = 1, e = . (15) 

The marginal operator is the local energy density, which 
keeps its anomalous dimension x™^ = 1, independently of 
the value of J or ft,. We note that nonuniversal interface 
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critical behavior at a defect plane can be found in the 
three-dimensional n- vector model in the limit n ~^ oo, 
which has been explicitly calculated,— 

D. Two different subsystems 

If the two subsystems have different four-spin couplings 
e^"' < e^^\ one can no longer define a bulk system fixed 
point. However, the ordinary and extraordinary inter- 
face fixed points still exist. The stability analysis of the 
ordinary interface fixed point can be performed along 
the lines of Sec lIII CT| leading to an interface exponent 
i/i — l~Xm ~^Xm \ The ladder perturbation is irrelevant, 
i.e., the ordinary interface fixed point is stable (unstable) 
for e(~)-|-e(+) > (< 0). Through duality, as described in 
Sec lIII C"3| the same type of stability is expected to hold 
for the extraordinary interface fixed point, too. Conse- 
quently, the directions of the RG flows are analogous to 
the case of identical subsystems in Fig. [TJ just the role 
of the bulk fixed point is taken over by a new special in- 
terface fixed point, located at Jg = 0(1), which controls 
a special transition. The expected RG phase diagram is 
given in Fig. [21 

The stability or instability of the special interface fixed 
point requires that the scaling dimension of the local 
energy-density operator satisfies 

4"* > 1 for e(-) + e(+) < , 
x"^* < 1 for e(~) + e(+) > (16) 

at this fixed point. 

In the borderline case, e^"-' + e^"*"-* = 0, the perturba- 
tion is marginal at the ordinary and extraordinary fixed 
points. It is interesting to determine whether the inter- 
face remains marginal for any value of J, as it happens 
in the symmetric case. In the truly marginal case (i) the 
local magnetization exponent is a continuous function of 
the coupling: xj^* = x'^\J) [as in Eqs. ^ and P^ ] 
and (ii) the scaling dimension of the local energy-density 
operator has to remain constant: x"'* — 1. 

IV. NUMERICAL STUDY 

The calculation of the scaling dimensions associated 
with the interface, xj^' and xj."*, is based on a finite-size 
scaling analysis of the matrix elements of the correspond- 
ing operators, as indicated in Eq. The ground state 
of the system, with a length L for the two subsystems up 
to 86 (38) for the magnetization (energy) exponent, has 
been determined using the DMRG methodi^i In order to 
obtain a good accuracy, we have generally kept around 
TO = 150 states of the density matrix. 

The magnetization density has been calculated using 
symmetry-breaking boundary conditions, with the two 
types of spins held fixed at both ends, ctj.^ = r|.^ = +1. 
The magnetization density is determined on both sides 
of the interface when the system is asymmetric. 
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FIG. 3: Interface critical behavior between two identical 
Ashkin- Teller models with bulk four-spin coupling e = 0.5 
(relevant case). The upper part gives the magnetization ex- 
ponents and the lower part the energy-density exponents de- 
duced either from efnt on one side of the interface (left) or 
from efnt on the junction itself (right). 



For the energy density, we eliminate the regular contri- 
bution to the ground-state expectation value in Eq. ([9]) 
by taking the difference of the values obtained for the 
systems with free and fixed boundary conditions. Since 
the sign of the singular part generally changes when the 
boundary conditions are changed, a good precision can 
be obtained in this way;^^ As indicated in Eq. we 
calculate the energy density on the junction itself by tak- 
ing the ground-state expectation value (0|tT^(— l)cr^(l)|0) 
and on both sides of the interface, with (0|cr'^(±l)|0). 

From the values of the singular part of the matrix ele- 
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FIG. 4: As in Fig. [3]for the irrelevant case e = —0.5. 



FIG. 5: As in Fig. |3]for the marginal Ising limit e = 0. 



ment, say, mint(i) ^ L ^™ , at two different sizes, L and 
bL, we deduce effective exponents through two-point fits: 



lnmint(&L) - Inmint(i) _ int 



In 5 



(17) 



In order to obtain the same numerical accuracy for the 
different points, we keep the ratio b between neighboring 
sizes approximately constant. The effective exponents 
evolve toward their exact values when the mean size as- 
sociated with the two-point fit, {L) = L{b + l)/2, tends 
to infinity. 



A. Two identical subsystems 

We first check the validity of the phase diagrams given 
in Fig. [1] for the interface between identical critical sub- 
systems, the ladder defect in an otherwise homogeneous 
system. We have studied three values of the bulk four- 
spin coupling, e — 0.5, — 0.5, and 0, and calculated the 
interface magnetization and energy exponents for differ- 
ent values of the interface coupling J. The results are 
shown in Figs. 3-5. 

When e = 0.5 (Fig. [3]) the perturbation is relevant and 
the bulk fixed point unstable. For small values of the 
interface coupling, the flow is toward a free surface be- 
havior. For J — 0.2, the effective exponents tend to their 



surface values, x. 



int 



2/3 and either x] 



int 



when the energy operator is the surface energy operator 
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of one subsystem (ef^j) or a;™' = 2xf^ = 4/3 when the en- 
ergy operator involves the surface magnetization opera- 
tors of the two subsystems (ef^^). The effective exponents 
converge slowly to a;™* = and x™^ ~ 2, characteristic 
of an ordered interface, for the highest values of J. The 
interface exponents take the bulk values, x™^ =1/8 and 
Xg"* = 3/2, for an intermediate value of J, between 1.25 
and 1.5, where the flow is toward the (unstable in the J 
direction) bulk fixed point. 

For £ = —0.5 the bulk fixed point is stable and the 
effective exponents in Fig. 2] approach the bulk values 
xj^* = Xm = and a;™* = Xg = 3/2, independently of 
the value of the interface coupling. 

The Ising limit e = in Fig. [5] is a truly marginal 
situation. As expected, the interface magnetization ex- 
ponent is continuously varying with J. The extrapolated 
values are in agreement with the exact results given in 
Eq. (Ull): xj^*(J = 0.2) = 0.382, xj^*(J = 0.5) = 0.248, 
and a;J^*(J = 1.5) = 0.07, respectively. The interface 
energy exponent takes the bulk value xj."' — 1, which is 
necessary for a true marginal behavior at the line defect. 



B. Two different subsystems 

For an interface between two different subsystems, we 
start with the case where e*^"^ -|- 6*^+^ > 0, which corre- 
sponds to the RG flow in the upper part of Fig. [51 

The results obtained for the magnetization (energy) 
density exponents when e'^-' = —0.25 and 6^+) = 0.5 are 
presented in the upper (lower) part of Fig. [SI In accor- 
dance with the RG phase diagram, for small J {J = 0.2 
and 0.5), the effective interface exponents slowly ap- 
proach the surface magnetization exponent of the right 
subsystem xf^ = 2/3, whereas in the other limit (J — 1.5) 
they seem to converge to zero. According to the numer- 
ical results, the special interface transition takes place 
at « 1 where the magnetization exponent is close to 
Xm = 1/8. 

The energy-density exponents shown in the lower part 
of Fig. [Sj are greater for small and large values of J than 
for J = 1.0 ~ Js- This behavior is expected since xl = 2 
at the ordinary and extraordinary transitions, whereas 
the instability of the special interface fixed point requires 
x'f < 1. 

When e(-) -he(+) < 0, we are in the situation sketched 
in the lower part of Fig. [5| which was tested for e'^"-' — 
—0.5 and e(+) = 0.25. The numerical results are pre- 
sented in Fig. \7\ 

Here, too, the crossover effects are quite strong for the 
effective magnetic exponents shown in the upper part 
of the figure. For a small coupling J = 0.2, the effec- 
tive exponents remain close to the surface magnetization 
exponent of the e^^^ = —0.5 model, xf^ = 1/3, with 
a tendency to decrease at the largest sizes. The value 
of the coupling J, at the special interface fixed point is 
slightly higher than 1, where the effective magnetization 
exponents have the smallest finite-size corrections and 
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FIG. 6: Interface critical behavior between two Ashkin- Teller 
models with four-spin couplings e'~^ — —0.25 and e'^-* = 
0.5 (relevant case). The upper part gives the magnetization 
exponents calculated either on the left (full symbols) or on the 
right (open symbols) of the interface. The lower part gives 
the energy-density exponents deduced from ef^t on the two 
sides of the interface (left) or from ef^t on the junction itself 
(right). 



the extrapolated value is a little below Xm = 1/8. 

The stability of the fixed point of the special interface 
transition is related to the value of the energy-density 
exponent xj."*, which is shown in the lower part of Fig. [T] 
Except for J = 0.2, the effective exponents extrapolates 
to values larger than 1, in agreement with the stability 



„int 



analysis in Eq. p6)) . For J = 1.0, one obtains 
1.10(2), which corresponds to a crossover exponent j/j = 
—0.10(2). This small (negative) value of the crossover 
exponent explains the slow convergence of the effective 
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FIG. 7: As in Fig. [6] for e*"' = -0.5 and = 0.25 (irrele- 
vant case). 



FIG. 8: As in Fig. [6] for the marginal case with e^- ' = —0.25 
and = 0.25. 



DISCUSSION 



magnetization exponents in the upper part of Fig. [T] 

For the marginal situation where t^^'^ + = 0, we 
considered two cases, e^+^ = 0.25 and 0.5. The results 
are shown in Figs. [8] and [9l respectively. 

The magnetic exponents seem to vary continuously 
with J, without evidence of crossover effects at large 
sizes. The possibility that this system is truly marginal 
is supported by the behavior of the effective energy ex- 
ponents which, whatever the value of J, extrapolate to a 
value compatible with xj."* = 1. 



One interesting feature of the interface critical behav- 
ior in the AT model is that the local critical exponents 
are continuously varying with the strength of the junction 
when the sum of the four-spin couplings vanishes, even in 
the asymmetric case. Here, we discuss the possible origin 
of this truly marginal behavior. 

We consider a somewhat different setting, where the 
system is semi-infinite and consists of two subsystems 
with the shape of corners, —oo < z < Q,Q < y < oo and 
0<2;<oo,0<?/<oo, connected by a chain junction 
along the sides at z = 0. We are interested in the behav- 
ior of the generalized corner exponent a;^, measured at 
y = z = Q. Under a logarithmic conformal mapping, the 



9 



0.5 
0.4 
0.3 
K~ 0.2 
0.1 
0.0 

( 

1.4 
1.2 
.5 1.0 
0.8 
0.6 



(a) e*"=-0.5 e*^=0.5 o j=o.2 




0.0 



0.1 



0.2 

<L>" 



0.3 



0.4 



(b) 













(c) 7=0.2 
□ J=0.5 
oj=1.0 
A J=1.5 




0.0 0.1 0.2 0.3 0.4 0.0 0.1 0.2 0.3 

<L>~' <L>~' 



0.4 



FIG. 9: As in Fig. [6] for the marginal case with ' = —0.5 
and e^+' = 0.5. 



critical semi-infinite system is transformed into a strip 
with open boundaries at i = ±L/2 and a chain junction 
at i = 0. In the Hamiltonian limit, the strip Hamiltonian 
Tistr involves a transverse field /i at i = 0, whereas the 
two ends of the chain are free. In the following, we cal- 
culate the first gap AE{h) of TCgti, perturbatively for a 
small transverse field, and deduce the local scaling dimen- 
sion through the gap-exponent relation of Eq. PU)) . 
where 27r has to be replaced by tt, the actual angle in 
the mapping of the semi-infinite system. To calculate 
the gap, we first perform the duality transformation in 
Eq. (21) . The transformed chain has fixed boundary spins 
a.t i ~ ±L/2 and a (weak) defect coupling of strength 
h at i = 0. The first gap is given by the difference of 
the ground-state energies with antiparallel and parallel 



boundary conditions: AE{h) = E^^ — Eq . Actually, an- 
tiparallel boundary conditions are applied to one type of 
spin variables, say, a, while parallel boundary conditions 
are always applied to the r spin variables. To leading or- 
der, only the a spin variables contribute to the difference 
of the ground-state energies and we obtain 



where 



AE{h) = 2m)r>{L/2)m\+>{L/2)h 
= 2a(~)aW(L/2)- 

^^'^\l/2) ^ a(±)(L/2)-^-'^' 



h, (18) 



is the surface mag- 
netization in a critical chain with length L/2, when the 
spin of the same type on the other end is fixed in the up 
state. The second-order term vanishes since even contri- 



butions to i?Q^ and to Eq^ are exactly the same. Only 
odd powers of h are present. 

The leading behavior of the gap in Eq. p8)) depends 
on the value of Xi — Xm + Xm^ . For Xi > 1, the first 
gap vanishes faster than 1/L; thus, x^ = and the junc- 



tion is ordered. This happens for 



e(+) > and 



corresponds to the upper part of Fig. [21 On the con- 
trary, for Xi < 1, the gap has a decay slower than 1/L; 
thus, according to Eq. pHl) . the interface exponent x^* is 
formally divergent to leading order of the perturbational 
calculation. This indicates that the extraordinary inter- 
face fixed point with ft, = is unstable, a situation which 
corresponds to the lower part of Fig. |2l In the marginal 
case = 1, up to first order in h, the local exponent has 
the variation 



Ah + 0{h^ 

TTV.. 



(19) 



where the coefficients a^^^ are 0(1). 

In the Ising limit e'^-' — 0, 0*^=*=^ = 1 and, with 
the parametrization chosen for the quantum Hamilto- 
nian, Vs = 2. Thus we obtain x^^(h) = ^ + 0(h^), 
which is the leading contribution to the exact result 
a;Jj^(/i) ;= 1 - 1 arctan(/i~^). 



In the asymmetric marginal case, e 



(-) 



0, we 



also have a continuous variation of the leading contribu- 
tion to with h. We also expect a truly marginal local 
critical behavior in this case. This assumption is sup- 
ported by the fact that the second-order term of the ex- 
pansion is vanishing due to symmetry. In the marginally 
relevant or irrelevant cases, the second-order term of the 
expansion is usually diverging as log however, for 
a marginally irrelevant perturbation, the singular terms 
are expected to sum up to a regular contribution. In our 
case, we expect a truly marginal behavior and continu- 
ously varying local scaling exponents also for the ladder 
junction studied numerically in Sec. II VI 

There are other systems from which similar compos- 
ite critical systems can be built and for which a truly 
marginal interface critical behavior could be obtained. 
Let us mention the 2D XY model with different tem- 
peratures, say, T^"^ and T^^\ both lower than the 
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Kosterlitz-Thouless temperature. Another example is 
the XX Z chain with different anisotropies on the two 
sides of the junction. One may notice that the AT 
Hamihonian can be transformed into a staggered XX Z 
model through a duality transformation of the r spins 
followed by a duality transformation on all the spinsi^ 
Finally, let us mention the Potts model in the Fortuin- 
Kasteleyn representation, for which the number of states 
q becomes a continuous parameter. Taking two subsys- 
tems with q states on one side and 4 — g states on the 
other, since xf^{q) = 1 — (2/7r) arccos(y/^/2),— one has 
^mil) + ^m(4 ~ = 1- It follows that the junction is 
a marginal perturbation at small coupling. In this case, 
however, the central charges of the conformal field theory 
at the two sides of the junction are different; therefore. 



it needs further investigations to decide if the perturba- 
tion remains truly marginal for any coupling strength, as 
observed in the symmetric Ising limit q = 2. 
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